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Abstract We consider quantum analogs of the relativistic Toda lattices and 
give new 2x2 L-operators for these models. Making use of the variable separation 
the spectral problem for the quantum integrals of motion is reduced to solving 
one-dimensional separation equations. 



1. Introduction 

The relativistic Toda lattices (RTL's) were originally introduced in [9] where the 
classical nonperiodic case has been solved by means of an explicit action-angle 
transformation with the help of the N xN matrix Lax representation for the model. 
For this RTL, which is naturally associated to the root system Ajv_i, the Lax 
representation in terms of N x N matrices appeared to have more complicated 
structure in comparison with the Lax pair for corresponding non-relativistic Toda 
lattice [8]. In [1,2] there was introduced a Lax triad for the model, three matrices 
of which have simplier three-diagonal form. With the help of this representation 
the direct and inverse problems for the system have been solved and the separation 
of variables has been carried out in full analogy to the classical work [4]. 

In recent works [12-14] there was proved an equivalence, due to canonical trans- 
formation, of the RTL's for general classical root systems and discrete time Toda 
lattices. There was given also an orbit interpretation of the Lax representations 
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introduced 2x2 L-operators which satisfy the Sklyanin quadratic algebra with 
non-standard trigonometric r-matrix. 

In the present paper we construct, for all the RTL's, new 2x2 quantum L- 
operators obeying the qudratic .R-matrix relations with the standard trigonometric 
.R-matrix. We give also the separation of variables for the A^-i RTL which is a 
generalization to relativistic case of the variable separation procedure given in [10] 
for non-relativistic Toda lattice. 

The structure of the paper is as follows: in the Section 2 we recall the known facts 
about the classical RTL, which will be needed in the sequel. In the Section 3 there 
is introduced new 2x2 L-operator for the A^-i RTL, which is a nonsymmetric 
L-operator of the lattice sh-Gordon model, and there is given also the variable 
separation procedure indicating that our separation variables are equivalent to the 
ones introduced before [1]. In the Section 4 we consider the quantum An-i RTL 
and construct the separation variables and corresponding separation equations for 
it. The final Section 5 describes the construction of the monodromy matrix for the 
quantum RTL's for general classical root systems (including affine cases). 

2. Description of the model 

The relativistic generalizations of the classical Galilei-invariant Calogero-Moser sys- 
tems and Toda lattices are characterized by the time and space translation gener- 
ators 

T=\(S+ P=\{S+ -S-), 

(1) " N 

where S± = ^ e ±dj Vj {qi, . . . , q N ) , 
and the boost generator 



N 



(2) B = -J2qj- 

3 = 1 

In this equations 9j and qj are canonically conjugated momentum and coordinate 
of the j-th particle. One obtains a representation of the Lie algebra of the Poincare 
group if and only if Vj satisfy the functional equations [9] 



VidiVj +V j d j V i =0, i^j 

d_ 

dqj 



J = l " ^ 

There are known two solutions to this equation [9] : 

(4) l.V j (q 1 ,...,q N ) = l[f(q j -q i ), f =a + bp(q) 

where p is the Weierstrass function. 

(5) 2. Vj (qi , . . . , q N ) = / (qj-i - qj) f (q~j ~ qj+i ) , 

P — /I , „2 „Q \l/2 „ r- TD 
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The first case is called relativistic Calogero-Moser system and the second case is 
referred to as RTL. 

As phase space we have 

N 

(6) O = {(9 , q) G R 2N } , u = d9 J A d Qj ■ 

3 = 1 

We will consider RTL as the integrable system of N one-dimentional particles 
with the Hamiltonian 

^ r i 1/2 

( 7 ) H = 2^ ex P ( d J + ex Pfe _ ) ] [ 1 + exp(g i+ i -qj)] 

3 = 1 

The RTL is open (non-periodic) when M = N — 1 and periodic in the case of 

(8) M = N, qj+kN = Qj > 9j+kN = where j = 1, . . . , AT and k e Z . 

In this section we consider only periodic RTL. 

According to [12,13] there is an equivalence between the RTL and the discrete 
time Toda lattice. It appeared to be more natural to introduce new canonical 
coordinates (pj ,qj) connected to the old coordinates (6j , qj) as follows: 

(9) Pj =6 J + 1 -ln( 1 + eXP ^-^- 1 
{) 3 3 2 \l + exp(q J+1 

Transformation (9j , qj ) — > (j>j , qj ) is a canonical transformation. 
In new variables the Hamiltonian (7) has the form 

N 

(10) H = Y1 ex Pfe ) 1 + ex P(<?j+i - Qj ) ] • 

One can introduce a set of new variables Cj , dj and /j [1,2] 

(11) cj = exp(pj - q~j + q j+1 ) , dj = ex.ppj , /? = q . 

Then the equations of motion given by the Hamiltonian (10) have the following 
form 

( 12 ^ d 3 = - d 3 ( C 3~l ~ C j) > 

dj = c 3 (d j+1 - dj + cj+i - Cj-_i) , 

and can be rewritten as a compatibility condition for the two linear problems: 

(k 2 - dj) (pj + k (fj cp j+1 + fj- x (pj-! ) = , 
(!3) -1 2 1 

„.l At IX A \t ■ „„ „„11„J T3„U Auu;„„ f, -t; 
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Let us give the equations (13) in the matrix form introducing two NxN matrices 
C (k) and M (k) 



(14) 



C(k)(f) = 0, j>=-M(k). 



The equations of motion are thereby written in terms of the Lax triad (£, M., J-) 
[1,2]: 



(15) 



d -c 

dt 



C,M\+F£, 



where the symmetric matrix C has the form 



/k 2 -di kft 



(16) 



kfi 
V kf N 







k 2 -d 2 kf 2 



kf N \ 




kf N -i 
■ ■■ k /at-i k 2 -d N / 



and the matrix M. looks like 



M 



1 

2 



/Ci - Cat + /c 2 
-kfi 



kfi 

C 2 -C! + k 2 





fe/2 



-kfN 




The matrix JF is diagonal and has zero trace: 



(fl-f 2 



N 



fi fi 

f 2 - f 2 J 

JN JN-1 / 

From the Lax representation (15) there follows that 



kfN-i 

-fc/jv-i cn - cjv-i + k 2 J 



\ 



(17) tr£ £ = tr^"=0, 

and, therefore, using the Abel identity 



(18) 



dlndetX x ^ 

— tr JC 



dt 



dt 



one concludes that the det (£) is the generating function of the integrals of motion 
for the RTL. The Lax triad (15) and the related questions were studied in [1,2]. 
We are giving below one Proposition from there. 

Proposition 1. 

a). The determinant and all principal minors of the matrix C (k) are polynomials 

u2 



b). The zeros of determinant and of all principal minors of the matrix C{k) are 
simple and positive real. 

According to [3,7] it is possible to rewrite the equations (13) in the matrix form 
with the help of the 2x2 matrices Lj (k) [12] 



(19) 



4-w = (/- exp(p - 



3) i 



-k exp((/j ) 




exp(pj- -qj), 

The matrix Cn (k) is related to the 2x2 monodromy matrix Tn (k) such that 

JV 



(20) det C N (k) - 2 = tr f N (k) = tr ] J Lj (k) = tr 

i=i 



An B n 
Cn D n 



Matrices Lj (A) and monodromy matrix Tjv (A) satisfy the Sklyanin quadratic al- 
gebra [3,10] 



(21) 



(T (1) (A), M } = [f (A/p) , f W (A) T( 2 ) (^ 



where (A) = T(A) <g> J and T^ 2 ) (//) = / <g> T(//) , with the following non- 
standard r-matrix 



(22) 



f{X/ti) 



- - M 2 



fa 








-1/2 





c 


Vo 







A /i 



\ 

c 
1/2 



rewritten in the form 
(23) {T( A) ® T (p) } = [r (A/p) , T ( A) ® T (p) 



+ 



s,T(A) (g) T(/x) 



In equation (23) r(A//i) is the standard r-matrix of the XX Z type and s is the 
matrix of scalars: 



(24) 



fa \ 

c 

c 

VO a ) 



fO 0^ 

-1/2 

1/2 

Vo ) 



In the quantum case the equation (23) goes to the nonstandard quadratic relation 

(25) R x (A///) T« (A) (p) = (p) T« (A) R 2 . 

We can express the entries of the monodromy matrix T (A) in terms of the 
determinant and principal minors of the matrix C(k) in the same way as it was 
done in [7] for the non-relativistic Toda lattice. Then using the Proposition 1 we 
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Proposition 2. 

a) An (A) = A^ + . . . is a polynomial of degree N on the A = k 2 where: 

N 

l 7V (0) = (-l) Ar exp(P), P = Y J Po- 

i=i 

b) B N (A) = k ((3X N ~ 1 + ...), C N (A) = k (tA^" 1 + ...) l D N {\) = 5X N ~ X + .... 

c) Zeros of the polynomials A (A), B(X)/k, C(X)/k, and D(A) are simple and posi- 
tive real. 

These facts will be used in the next section for doing the separation of variables. 
We will give also a new 2x2 L-operator and a new monodromy matrix for the 
RTL which are closely connected to the previous ones, but satisfy the standard 
XX Z-tjpe quadratic r-matrix algebra. 

3. Classical mechanics 

It is well-known that the 2x2 L-operator for the non-relativistic Toda lattice can be 
obtained by contraction from the L-operator of the XXX-tjpe Heisenberg magnet. 
For the XXZ-type Heisenberg magnet one can use the analogous contraction and 
construct the L-operator of the form: 

(26) L J (u)=( 2sinh{u - P i ) - eXp( ^ 

\ exp(-^ ) 

where the pj and qj are the canonically conjugated momentum and coordinate of 
j-th particle. This L-operator obeys the fundamental Poisson brackets (21) with 
the standard r-matrix r{X/ ji) of the XXZ-type (24) (A = e u ) and it can be thought 
of as a nonsymmetric analog of the L-operator for the lattice sh-Gordon model [6]. 
The L-operator (26) and the monodromy matrix Tjv (u), 

N , x 

(27) T N (u) = f[L j (u)=(£ B D ) (u), 

j=i \ / 

have a simple connection to the Lj (fc)-operator (19) and the monodromy matrix 
T (k) (20) through the following change of variables: 

Qj =9j - y , k = exp(w) ; 

(28) Lj (u) = k~ x exp(— ^-) L (u) , 

T N (u) =k~ N exp(-^)f N (u), 

where P = XljLi Pj ^ s ^ e total momentum of the system which is an integral of 
motion. 

We can associate the following N x N matrix L (u) to the monodromy matrix 
T N (u) (cf. with [7]): 

/2sinh(u -pi/2) e 2 i ei N \ 

e 2 i 2sinh(-u — p 2 /2) e 32 
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where e jn = exp(^— - — -). The matrices L{u) and C(k) are connected by the 
formula 

(29) L(u) = UC(k)U 

where the matrix U has the form 

U = diag(exp(-|), exp(-^), . . . , exp(-^) ) . 

Due to the relation (29) one can construct the Lax triad for the matrix L (u): 

dL{u) 



dt 



L , M 



+ FL , 



where the matrices M and F are expressed in terms of the known matrices Ai and 
F(see (15)-(16)) 

M = U~ X MU + U~ 1 U , 
F = U- X MU -UMU' 1 +f+2U- 1 U . 

Notice that in the limit to non-relativistic case the Lax triad L , M , F turns into 
well-known Lax pair for the usual Toda lattice [4,8]. 

The separation of variables for the periodic RTL was originally given in [1]. 
Let us carry out this procedure within the r-matrix method with an intention to 
generalise all the results afterwards to the quantum case. 

We begin with the monodromy matrix T (k) (20)-(21), since due to the Propo- 
sition 2 the entries of it have simple dependence upon the spectral parameter k. 
Let us take as new variables \ij the zeros of the off-diagonal entry C (A) (for such 
chosen, the variables \ij coincide up to a factor with the ones introduced earlier [2] 
and are zeros of the Baker- Akhiezer function): 

(31) C(fjLj) = 0, j = l,...,N-l. 

By the Proposition 2 these variables are positive real and independent. Further, 
define a set of the variables as follows: 

(32) ur=A N ( tij ), v+=D N ( H ), j' = l,...,iV-l. 
The variables \ij , , j = 1, . . . , N — 1 obey the Poisson brackets: 

d detT(A) 



(33) {v, ,^} = Sj 



dX 

u+ vj = det f (fij ) . 



These brackets follow from the fundamental Poisson brackets (21) and can be cal- 
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Introduce two more variables 

/j N = exp(p N -q N ), Vn = ex P( ±P )- 

The set of the fij , up j = 1, . . . , N is complete, i.e. any dynamical variable can 
be expressed through them. To prove this let us restore the monodromy matrix 
T (k) in terms of the new variables: 



A N (\)=xu- n (a-^) + e n ^t v . 

j = j = % = ^ ^ J 



where ^0=0 and z/ Q =1, 
(34) 

N-l N-l . AT-1 

^v(A) = £ ±^-v}, (7 Ar (A)=^fc (A-/*,) 

J = l, = l / '' , = 

» ^ J 

(recall that A = k 2 ). The entry B (k) of the matrix T (k) can be obtained using 
the identity det f (k) = k N exp( -P) . 

Now we introduce another set of variables yj , Yp using the monodromy matrix 
T(«)(26)-(27): 

(35) C( % ) = 0, r- = ^( % ), r+ = J D( % ), j = l iV-1. 

The matrix entry C (u), considered as a function on the spectral parameter u, has 
infinite number of zeros, but using the connection between the matrices T (u) and 
T (k) (28) we can prove that these zeros have the following structure 

(36) yj n = yj + 2% in = In \ij + 2tt i n , n = 0, 1, 2, . . . , 

hence the variables yj , j = 1, . . . , N — 1 can be chosen to be simple and real. 
Variables yj , Yp obey the Poisson brackets: 

{Vj,Vn} = {Y?,Y±} = 0, YpYp=i, {y J -,y B ± } = ±tf in y B ± . 

One can restore the monodromy matrix T (u) in terms of these variables. 

The set y 3 , Yp gives the separation variables, since by their definition (35), they 
obey the equations 



(37) 



t r T N ( Vj ) = Y ■+ + Y - , j = 1 , . . . , N - 1 . 
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4. Separation of variables in the quantum case 

In this section we apply the general scheme [10] of the quantum separation of 
variables to the quantum An-i RTL model. 

Quantum L-operator for the non-relativistic Toda lattice [3,10] has the form 



Lj (u) 



u 

,exp(-^ ) 



— -exp(^ 







For the RTL the quantum L-operator and the monodromy matrix TV (u) are 



Pi 



L, (u) 



(38) 



2sinh(-u 
exp(-qj ) 



-exp(ty) 



n f A B \ 

T N («) = \\ Lj (u) = I c D j , where [pj ,q n ] = ir]5 jn 

7 = 1 ^ ' 



and obey the quadratic .R-matrix algebra 

(39) R(u-v) TW (u) T< 2) (v) = (v) (u) R(u-v), 

where (u) = T(u) <g> 7 and T^ 2 ) (v) = I ® T(u) and il(u) is the following 
trigonometric solution of the quantum Yang-Baxter equation [3,6,10]: 



(40) 



i2(u) = 



fa{u) \ 

b(u) c(u) 

c(u) b(u) 

V a(u) J 



where 



(41) a (u) = sinh(w + rf) , b (u) = sinh-u , c (u) = sinh?? . 

It is natural now from the above formulation of the problem to call the RTL a 
(/-deformation (q = e v ) of the ordinary Toda lattice. 

Separation of variables for the quantum RTL is taken in a complete analogy to 
the separation for the ordinary Toda lattice [10]. One considers the operator roots 
of the equation Cn (u) = 0. They are N — 1 Hermitian operators yj which give 
quantum analogs of the corresponding real simple zeros in the classical case (36). 
Let us introduce the operators Y^, j = 1, . . . , N — 1 by the rule: 

(42) Yf = A N (u => yj ) , Y+ = D N (u => y 3 ) . 

These operators are defined by the values of the functions A (u) and D (u) in u = yj. 
The substitution of operator values for u will be defined correctly if one prescribes 
some rule for operator ordering. Here the ordering of y's to the left chosen. We will 
call it substitution from the left and denote u => yj. The operators yj , Y^ , j = 
1, . . . , N — 1 have nice commutation relations 

L. „. 1 _ f^± v±I _n ^±i^=F _i L. v±1 _ ■ x v± 
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which follow from the basic commutation relation (39) [10]. 
Let we consider the spectral problem 

trT(u) ip = r (u)ip 

then substitute u => yj (from the left) and use the definition of Y ? ± . The resulting 
set of equations has the form 

r(y J )^(y J ) = (A.Y-+A + Y+)i;(y J ) 

where A± (u) give a factorization of the quantum determinant of the monodromy 
matrix [10] A(u) = det q T (u) = A + A_ = 1. Let A_ = i~ N , A + = i N . In 
^-representation the operators are acting as shift operators and the separation 
equations take the form 

(43) t ( yj ) ip ( yj ) = i N ifj( yj + irj) + i~ N *p (yj - it/) 

while the full eigenfunction can be represented in the factorised form (quantum 
separation of variables) 

N-l 

ip(yi,... ,y N -i) = IJ ^(Vk)- 
k=i 

The separation equations above are — 1 one-dimensional finite- difference multi- 
parameter spectral problems. 

We can rewrite equations (43) as 

N 

( In -2k exp( (N - 2k) y 3 ) ) ^ ( Vj ) = 

k=0 

= i N ip(yj + ir]) + i~ N ip(yj - irj) 

where I±k are the eigenvalues of the quantum integrals of motion. For ordinary 
Toda lattice these equations have the form 

N 

J2{hyj)ii(y j ) = i N ip(y j +iv) + ^ N ^(yj -*v)- 
fc=i 

As for a solution of the above separation equations, there are more questions till 
now than answers, although in the non-relativistic case one can use the Gutzwiller 
algorithm. The RTL case provides even more difficult problem, for example, for 
N = 2 the separation equation 

(cosh p + E ) ip = sinh x ip 



(44) 
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5. MONODROMY MATRICES OF THE QUANTUM RTL 
FOR GENERAL CLASSICAL ROOT SYSTEMS 

The monodromy matrices for such systems obey, as well as in the non-relativistic 
case [7,11], the reflection equations 

1 2 

R{u-v) U- (u)R(u + v- 7]) U- (v) = 



(45) 



2 1 

U- (v)R (u + v-7])U- (u) R(u-v), 

1 2 

R (-u + v) Ul 1 (u) R(-u-v-rj) (v) = 

2 1 

(v)R (-u-v- rj) U% (v) R (-u + v) , 

1 2 

where U± = U± <8> / \U± = I <8> U± and are matrix transpositions in the first 
and second spaces, respectively. 

We have the following rules for constructing the monodromy matrices T(u) here 
[11]: 

U-(u) = T_(u)K_ {u-^)Tl\-u), 

(46 ) U l («) = T l («) K^+j) ("«) 

T(u) = U + (u) U- (it), for the closed lattices (interval) 

T(u) = U- (it), for the open lattices (semi-axis) 

where the T± (u) are matrices obeying the equation (39) and the K± are simple 
solutions of the reflection equations (45). The generating function of the integrals 
of motion is [7,11] 

t(u) = trf(u) = tr U+U- 

= tr K + (u+^)T(u)K_ (u-^)T^(-u), 

M 

t («) = n Lj («) . 

For the root systems and Cat, M = N and matrices lf± are scalar solutions 
of the reflection equations (45); for the root system D N , M = N — 2 and matrices 
K± depend on the dynamical variables. 

We list below the matrices K± and corresponding Hamiltonians for the quan- 
tum generalized RTL's (we give also the corresponding K- matrices for the non- 
relativistic Toda lattices): 

Bn and series: 

for RTL's 



(47) 





+ 


ft 


sinh u 


cosh u 




7i 




ajv 




Pn 


sinh u 




cosh u 



-7i 

Ql Pi 

sinhu coshu 

In \ 

rv*7 . Bkt I 1 
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for non-relativistic Toda lattices 



K_ = 



-U 



K+ = 



Pi U Oil 

The Hamiltonian of the corresponding RTL's is 
(48) H = H +J 1 +J Nl 

where we use the notations 



«7V —/3nu 

—U CtJV 



N 



JV-1 



(49) H = 2 cosh( Pj ) + 2 Y 



exp(q J+1 -qj)- cosh(^— J -) 



(50) 



Ji=2exp(^i) «icosh(y) +/3isinh(y) + 71 exp(2<?i ) , 
J N = 2 exp(q N ) a N cosh( — )+(3 N smh( — ) + ~f N exp(-2q N ) . 



We have here RTL's associated to the following root systems: Cn («i = 
Pi =0), Bn {oti = Pn = Pi =0), for the open lattices and Cff («i = cxn 
0), A^ N (Pi = &n =0), D$ + i (Pi = Pn = 0) for the closed lattices. 

D N : 

Denote Fj = cosh (pj + qj ) and Gj = cosh (pj — qj ) then 

^ ( exp (u) Fi + exp (—u) Gi sinh 2 qi — sinh 2 (u) 

sinh 2 (u) — sinh 2 pi exp (u) Fi — exp (— u)Gi 



(51) 



K + = 



/ exp (u) Fn — exp (—u) G 



N 



sinh 2 (u) — sinh 2 pj 



sinh qi — sinh (u) exp (u) Fn + exp (— u)Gn 
For ordinary Toda lattices 



ucoshqi — pi sinh^i sinh 2 qi 

-ucoshgi + pi sinhgi 



u 2 - pj 



. u cosh q N +Pn sinh q N u 2 - p 2 N 

1 sinh 2 qN u cosh qN — pn sinhgjv J 



The Hamiltonian now is H = H + Ji + Jn where H is given by (49) and 
(52) 



Ji =2 exp(<?i + q 2 ) cosh(^-^-^ 2 -) + exp(2g 2 ) 



70/ ^ urP N +PN-1 \ 1 / n \ 
J N = 2 exp(-q N - qN-i ) cosh( ) + exp(-2gjv_i ) , 

We have here RTL's associated to the following root systems: , B$ (Pn = 

„„ j A ( 2 ) /„. _ n\ 
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